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The neutron index of refraction is general- 
ly derived theoretically in the Fermi 
approximation. However, the Fermi 
approximation neglects the effects of the 
binding of the nuclei of a material as well 
as multiple scattering. Calculations by 
Nowak introduced correction terms to the 
neutron index of refraction that are quad- 
ratic in the scattering length and of order 
1 0~ fin for hydrogen and deuterium. 
These correction terms produce a small 
shift in the final value for the coherent 
scattering length of H2 in a recent neutron 
interferometry experiment. 
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1. Neutron Optics Theory 

The basis of neutron optics is the Schrodinger equa- 
tion for the coherent interaction of the neutron with 
matter using an optical potential. It is well known that 
the wavefunction satisfies the Lippmann-Schwinger 
scattering equation [1] 



V/(r) = |/:) + Gv(r)v/(r) 



(1) 



where v^r) is the wavefunction, \k) is the incident 
wave, G is the one-body Green's function, and v is the 
optical potential. It is related to the /-matrix by 



/ = v(r) + /Gv(r). 



(2) 



Thus, given a form for the /-matrix, one can determine 
the optical potential and use this to solve the 
Schrodinger equation for the form of the wavefunction. 



The form for the /-matrix used in the Fermi approxi- 
mation, the standard approach to solving for the optical 
potential in the theory of neutron optics, is 



'=Io. 



where 



m i 



(3) 



(4) 



In the above equation, j denotes the atomic species, m 
is the neutron mass, bj is the coherent scattering length, 
5 is the delta function, r is the coordinate of the neu- 
tron, and R^ is the coordinate of the atomic scattering 
center. Using this in conjunction with the Bom approx- 
imation, / = V, the optical potential then becomes 



v(r) = ^XA^,Z>, 

m j 

where A^ is the atomic density 



(5) 
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This leads to the familiar form of the index of refrac- 
tion for the Fermi approximation, 



„-l-^XA^,. 



(6) 



However, this approximation disregards atomic binding 
and multiple scattering. Thus Eq. (3) must be modified 
to take these effects into account. The next order term 
in the multiple scattering expansion of the /-matrix is 

In general this correction term to the /-matrix is diffi- 
cult to evaluate. However, in the particular case of 
molecular hydrogen and deuterium gas, one can per- 
form this calculation in the rigid rotor approximation. 
The two correction terms for H2 gas are 



iN ( Inh 



(InyM m 






and 
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(9) 



due to the binding potential and multiple scattering 
respectively. These potential terms may be expressed in 
terms of the index of refraction, 

AkNE 

-— ^(i+4,+4.). (10) 



Table 1. Average corrections to the refractive index for H2 gas 



Temperature (K) 



A = A^, + Ah, 



100 
300 



2.0 X 10"- 
1.9 X 10"- 



Institute of Standards and Technology (NIST) Center 
for Neutron Research (NCNR) Interferometer and 
Optics Facility using a triple-blade, single perfect crys- 
tal silicon neutron interferometer. See Fig. 1 . A more 
detailed description of the method, experimental 
arrangement and procedure for the precise determina- 
tion of coherent scattering lengths using neutron inter- 
ferometry can be found in Schoen et al. [2]. 
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Fig. 1. A schematic view of the Si perfect crystal neutron interfer- 
ometer with gas cell. 



The phase shift that occurs between the two coherent 
neutron beams due to the presence of the H2 gas along 
one of the paths is given by 



If one assumes a rigid rotator approximation for the 
rotational excitations of the hydrogen molecule, the 
potential terms can be numerically integrated. Nowak 
performed these calculations [3], and his results are 
seen in Table 1 . 



2, Experimental Procedure 

The measurement of the real part of the n-p coherent 
scattering length was performed at the National 



A(l) = in-\)kd=-Nb^Xd, 



(11) 



where A0 is the phase shift, n is the real part of the 
index of refraction, k is the wavevector, A^ is the atom- 
ic density, X is the neutron wavelength, and d is the path 
length that the neutron beam traverses through the H2 
gas. From this relationship the coherent scattering 
length can be determined. 

A secondary sampling method is used to measure the 
phase shift. This is accomplished by positioning a 
quartz phase shifter across the two beams traversing the 
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interferometer and rotating the angle, 5, about the 
mounting axis. The intensities of the beams that arrive 
at the ^He detectors, labeled O beam and H beam as in 
the diagram, are a function of the phase shifter angle, 
and are given by 

I^{5) = A^+Bcos{Cf{S)+(^,) (12) 

4(5) = A^Bcos{Cf{5)+(^, +7t). (13) 

The function y( 5) is given by the path length difference 
between the two beams travelling through the phase 
shifter. Aq, A^^, B, C, and 0o are parameters that are used 
to fit the data. 

A gas cell with two chambers, specifically designed 
to minimize the phase shift due to the presence of the 
cell itself by spanning both beams, held one chamber 
filled with 99.999 % chemically pure H2 gas from 
Matheson Tri-Gas, and the other chamber was evacuat- 
ed. This cell was placed in the interferometer, after its 
mounting had been aligned using a quartz alignment 
flag, such that one beam encountered a vacuum along 
its path length, while the other encountered the H2 gas. 
The phase shift caused by the gas was determined by 
collecting interferograms. See Fig. 2. First an interfero- 
gram was measured with the gas cell in the interferom- 
eter, and the angle of the quartz phase shifter was var- 
ied. The phase was determined from the data fitting. 
This was then repeated with the gas cell out of the inter- 
ferometer, to obtain the interferometer offset phase. 
The difference between the two phases is due to the 
presence of the gas and the cell. The phase shift due to 
the cell itself was determined by holding both chambers 

m 
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Fig. 2. An example of an interferogram in which the intensity 
changes as the phase angle 5 is varied. Data is shown with both the 
gas cell in the interferometer and outside. 
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evacuated, and this value was subtracted from the gas 
and cell measurement to yield the phase shift due to the 
gas alone. This procedure of measuring the phase shift 
due to the gas was repeated approximately 350 times. 

The atomic density was determined using the virial 
equation up to the third virial coefficient. The second 
and third virial coefficients were obtained from the 
compilation by Dymond [4]. In addition to the correct 
coefficients, the virial equation also requires accurate 
knowledge of the temperature and pressure. The tem- 
perature was monitored and recorded using 1 00 Q. plat- 
inum resistance thermometers, and the pressure was 
measured using a silicon pressure transducer as in pre- 
vious scattering length measurements [2]. 

The thickness of the gas cell was measured to be 
(1.006 ± 0.0001) cm by the NIST Precision 
Engineering Division Coordinate Measuring Machine 
[5]. This thickness was corrected for thermal expansion 
and contraction due to temperature fluctuation (a = 2.5 
X 10"^ °C~^ for aluminum) for each individual run. The 
effect of the gas pressure (-12 bar) on the thickness of 
the cell was negligible. 

The wavelength of the neutrons traversing the inter- 
ferometer was measured using a pyrolytic graphite 
(PG 002) crystal. This analyzer crystal was placed in 
the H-beam of the interferometer and rotated about 
until both the symmetric and anti-symmetric Bragg 
reflections were found. Bragg's law, X = 2d sin6^, was 
then used to determine the actual wavelength. For this 
experiment, X = (0.2713050 ± 0.0000085) nm, and, in a 
separate test, the stability of the wavelength over time 
was shown to be 0.001 %. 

The value of the coherent scattering length was cal- 
culated for each data set on a run-by-run basis, and then 
averaged to obtain a final value of Z)h = (-3.7458 ± 
0.0020) fin. This result uncorrected, for the correction 
terms discussed above, differed from the world average 
of previous results, Z>h = (-3.7410 =b 0.0010) fin, by 
approximately 2.3 (7. The highest-precision previous 
measurements of the coherent n-p scattering length, 
conducted using a gravity refiactometer using neutron 
reflection fiom liquid hydrocarbons, require a negligi- 
ble correction for higher order contribution to the t- 
matrix due to the larger mass of the molecules. From 
Eq. (10), it is observed that the Nowak correction term 
simply adds to the nuclear scattering length. Therefore, 
in determining the n-p scattering length, the observed 
scattering length for the neutron-hydrogen system was 
simply divided by (1 +z\^^ + z\^J. Our result became 
^«/; = (-3.7384 ±0.002) fin, assuming no error in the 
theoretical calculation. This agrees with the previous 
experimental values (see Fig. 3). A higher precision 



261 



Volume 110, Number 3, May- June 2005 

Journal of Research of the National Institute of Standards and Technology 



-J.M 


I ' 1 1 r ' M 

V 


'I'll 


- 




n lotal Reflection 






- 


-3.fiR 


o TransTmssion 
■^ Gravity Rcflcctoinctiy 
V TnterferoTTietry 
Current Weighted Average 




- 


1 --^^72 


I&t Standard Uncertainty 








O 1 


™„.*^«:_.-w,_;ff — 


... 


-3.76 




: 1 




- 


-3.80 


■ r 

1 


< . 1 , 1 


1 - J 


■' 



■ \.m J 



1950 1960 1970 1980 1990 20[)0 

Year of MeasureiTient 

Fig. 3. A plot of the bound coherent neutron scattering lengths for 
the n-p system along with reported uncertainties. Our result for H2, 
neglecting the Nowak correction term, is 2.3 a away from the world 
average of experimental results. Inclusion of the theoretical correc- 
tion brings our value into closer agreement, as can be seen. 



measurement of the n-p scattering length using H2 gas 
at different temperatures and ortho-para ratios could be 
used to experimentally isolate these correction terms. 
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